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Abstract. The thermodynamics of massless ideal gas of overlap quarks has been 
investigated numerically for both zero and nonzero baryon chemical potential /i. While 
the parameter M has been shown to be irrelevant in the continuum limit, it is shown 
numerically that the continuum limit can be reached with relatively coarser lattices 
for certain ranges of M . Numerical limitation of the existing method of introduction 
of chemical potential in the overlap formalism is discussed. We have also studied the 
energy density of free domain wall fermions in the absence of /i and estimated the 
extent of lattice in the fifth dimension L5 for which the overlap results are recovered. 
Interestingly, this value of L5 is also minimum for the same range of AI found in the 
overlap case. 
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1. Introduction 

The phase diagram of Quantum Chromodynamics(QCD) has been an important subject 
of study in the recent years and lattice gauge theory has emerged as a major tool for 
studying it. It is believed that spontaneously broken chiral symmetry is restored at high 
temperatures so it is very important to study the QCD thermodynamics with fermions 
having exact chiral symmetry on lattice. The most popular among such fermions are 
the overlap and the domain wall fermions. The overlap operator is highly non-local, 
involving inversion and square root of Dirac matrices, making it computationally very 
expensive to simulate. The domain wall fermions on the other hand are necessarily 
defined on a 5-D lattice so their computational cost is more than the standard 4-D 
fermions as the Wilson and staggerred, but is still less than the overlap. For the 
above mentioned practical reasons, we investigate numerically whether the irrelevant 
parameter M in such fermion operators can be tuned optimally to recover the continuum 
values of various thermodynamic quantities for free fermions with the smallest possible 
lattice sizes. This would give us an estimate of the optimum lattice parameters for full 
QCD computations with such fermions. In section 2, we compute the energy density of 
overlap and domain wall fermions in absence of chemical potential. It was shown [1] that 
the domain wall operator reduces to the overlap when the extent of lattice in the fifth 
dimension, L5 — >• 00. We try to find out at what L5, the overlap results are obtained. 
In the subsequent sections the energy density in the presence of chemical potential /i 
and the susceptibility expressions are computed and their approach to the continuum 
for different values of M are studied. 

2. Energy density in absence of chemical potential 

2.1. Overlap fermions 

The overlap Dirac operator [2] for massless quarks defined on a A^^ x A^^- lattice with 
spacing a and 04 in the spatial and temporal directions respectively, is given by, 

Dov = l + 'j5sgn{'y5Dw) , (1) 
where sgn denotes the sign function and 

T^ / \ ji^r '^rr 1+74 l~74r i 

Dw{x,y) = - Af 6.^^y - —[^x-i,y^— + -^—^x+ij 

- h^x-ly^^ + ^-^Uly] + [3 + ^]Sx,y. (2) 

The sign function of the matrix 75-DvK can be written as the sign of its eigenvalues [3]. 
These eigenvalues are positive and greater than zero in this case so the sign function is 
well-defined. The overlap operator hence can be diagonalised in the momentum space 
in terms of the variables 

3 

hi = — sin api , /i4 = sin 04^4 , /15 = M— ^(1— cosapi) (1— cos 04^4) 

O4 j=l '^4 
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Figure 1. (a)The variation of e/esB with M for fixed C = 4 for overlap fermions and 
(b)The domain wall fermion results matches with the overlap for > 14 with C = 4 
and M = 1.55. 



The energy density can be obtained as 

d\xiZ{V,T) 



V dT 



1 (dlndetDr, 



V 

,2/ 



2 ^ h\l - cos ap4) - h^hi sin ap4 i „ , \ / a\ 



where a = has been chosen after performing the a4-derivative. The same can be 
computed numerically keeping the aspect ratio C, = N/Nt fixed. The lattice expression 
so obtained was fitted to the ansatz A + B/N^ for a range of values of M. The higher 
order Nt dependent correction terms in the energy density expression due to finite 
size effects manifest themselves as an effective coefficient B{M). We seek those values 
of M for which B{M) approaches Ttt^/GO, its value in the continuum limit. The zero 
temperature part A, was subtracted from the energy density. The resultant e divided 
by continuum value esB was plotted as a function of Nt- From figure 1(a) we observe 
that the M-dependence is quite pronounced and for 1.55 < M < 1.60 the continuum 
limit is reached for the smallest possible lattice size(A^T' ~ 12). For the oft-favoured 
choice of M = 1, the convergence to the continuum is seen to be very slow. The energy 
density can be evaluated analytically by converting the sum over Matsubara frequencies 
ap4 = uj = (2m + Vjir/Nrp where m are integers, to contour integrals[l]. It has been 
shown that no M-dependence of the energy density exist in the continuum [?] but for 
observing that we need lattice sizes much greater than Nf = 32. 

2.2. Domain wall fermions 

The domain wall fermions [5j are defined on a 5-D Euclidian lattice in presence of a 
domain wall of height M. Taking the lattice spacing in the fifth direction ^ and 
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the number of sites Nr^ ^ oo for fixed L5 = N^a^, the effective domain wall operator is 

Ddw = (1 + + (1 - ^)7^ tanh(|f D^) . (5) 
Here is the bare quark mass and L5 is in lattice units. In the momentum space, using 
hi, h2, h^, /i4 and as given in ([3]), the energy density for massless quarks (m^ = 0) is 
given by 



ea 



4 



E 



4sinh[^]((-/i4/i5tt + h^l) cosh[^] + {hih^a + {hi + 5^)7 + 2h^s^p)) 



where 



cosh(^) - 2s(/i|p + /i57 + {h^p + h^a + 2h^{h^p + 7)) cosh[sL5]) sinh[^]) 

—2h^s sinh[2sL5]) 

a = -hi ,7 = 1- cos((2m + 1)h/Nt) , (7) 
h'' = hl + hl + hl + hl, s"" = h'' + hl , (8) 

(-sin^ap4 + /i5(l - cosap4))(-tanh^ + ^sec/i^^) 
^ " s^tanh^ ■ 

Following the same fitting procedure mentioned in the previous section and choos- 
ing M = 1.55 the Nt dependent part of energy density e was plotted for various values 
of L5. Clearly from figure 1(b) we observe that for L5 > 14, the energy density converges 
to the continuum case at Nt > 12. L5 ~ 14 is thus the optimum lattice size in the fifth 
dimension to converge to the overlap results. 



(6) 



3. Non-zero chemical potential 

In the presence of chemical potential fi = fia there is a critical value of fic beyond 
which the sign function is not always definedjl]. Hence we work with fi < fic, where the 
argument of sign function is positive definite and the calculation follows analogously as 
in section 2.1. The chemical potential can be incorporated as[6] and e~'^ multiplying 
1 ±74 terms respectively in ([2]) so as to cancel fi'^/a^ terms on the lattice [7j. The energy 
density expression on the lattice has the same form as in (jl]) but with r = /i/T = fiNx 
kept constant and h4 and h^ changed to[3]: 

CL Q 

h^ = sin (04^4— i/i) , /i5 = M— ^(1— cosapj) (1— cos(a4P4— z/t)) . (10) 

The energy density in presence of fi has been analytically shown to have the correct 
continuum limit The two observables of interest here are the change in the energy 
density, Ae(/i, T) = e(/i, T) — e(0, T) and the quark number susceptibility at /t = which 
in the continuum limit, are given by 

-^ = ^2+^ & X{0,T)=XSB = ^ (11) 
As shown in figure 2, both these quantities also approach the continuum limit faster for 
the same optimum range of M. 
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Figure 2. (a) Comparison of tlie lattice Ae/T^ with that for the continuum for r 
(b)The variation of x{0)/xsb vs Nt for different M and ( = 5. 



0.5. 



4. Conclusion 



The range 1.55 < M < 1.60 is best suited for the computation of different 
thermodynamic quantities of free overlap and domain wall fermions on the lattice as the 
continuum results are obtained on the smallest possible lattice size. This range seem 
to be quite universal even in the presence of fi. The L5 needed to obtain the overlap 
results starting from the domain wall formalism is also minimum for this range of M. It 
is anticipated that this range would remain roughly same even in the presence of gauge 
fields and thus this work justifies the use of M > 1 for faster computations even in full 
QCD with such fermions. 
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